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Abstract 



We investigate the dynamics of electrons in the vicinity of the Anderson 
transition in d = 3 dimensions. Using the exact eigenstates from a numerical 
diagonalization, a number of quantities related to the critical behavior of the 
^ ' diffusion function are obtained. The relation ri = d — between the corre- 

\^ • lation dimension D2 of the multifractal eigenstates and the exponent rj which 

■ enters into correlation functions is verified. Numerically, we have rj ~ 1.3. 

Implications of critical dynamics for experiments are predicted. We investi- 



\G \ gate the long-time behavior of the motion of a wave packet. Furthermore, 

' electron-electron and electron-phonon scattering rates are calculated. For the 

' latter, we predict a change of the temperature dependence for low T due to rj. 

^ ■ The electron-electron scattering rate is found to be linear in T and depends 

. on the dimensionless conductance at the critical point. 

a " 
o 

^ ■ I. INTRODUCTION 
> ' 

^ . The dynamics of non-interacting electrons in the vicinity or directly at the critical point 

^ ! of a metal-insulator transition is still an unsolved problem. In the last years, a number of 
works [p5|,p!|-p|j5|,p0|,|^,[r| revealed strong amplitude fluctuations of the eigenstates near the 



critical energy E^. The wave functions turned out to be multifractal objects, described by 
a set of generalized fractal dimensions Dg. On the other hand, investigation of correlation 
functions in combination with scaling arguments showed that on certain length and time 
scales the dynamics near Ec is governed by anomalous diffusion, described by an exponent 
rj. One possible definition for r] is the algebraic decay [|1^ of the static two-particle 



correlation function S{r, E,uj 0) near the critical energy E ^ Ec 



r 



S{r,E,u^O)^l^-^^j , r«e(i?) (1) 
where C,{E) is the localization length. In real space, S is defined as 

S{r,E,u;) := (y^6{E- - Ec,)6{E+ - Ep)^EAO)n-ir)^E+{r)nAO)) , (2) 

\a,/3 / 
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where E"^ := E ± hu/2, a,P label the eigenstates, and the brackets ( ) denote an impurity 
average. The correlation function S is connected to the diffusion function D{q,uj) via 

where in turn D{q,uj) is the generalization of the diffusion constant D in the metallic case 
and p{E) is the density of states. Indeed, Eq. (0) is a general expression for the two-particle 



correlation function compatible with particle conservation |TT|. The function D(q, tu) appears 
in the diffusion propagator 

P{q,uj):= . j. (4) 

—lu + D{q, uj)q^ 

which is the Fourier transform of the probability distribution P{r,t) describing the motion 
of an electron wave packet in a disordered system for t > 0, initially located at r = at 
time t = 0. 

In this work, we determine the correlation dimension D2 of the multifractal eigenstates 
and the exponent rj for the Anderson transition in 3 dimensions. The relation 

D2 = d-r], (5) 

which holds in d dimensions, is derived from the respective definitions of the exponents and 
verified numerically. We independently determine D2 via the box counting method applied 
to the spatial amplitude fluctuations of the critical eigenstates on the one side, and rj from 
a correlation function in energy space using a scaling form of D{q,uj) on the other side. 
Furthermore, the relation 5 = D2 = D2/ d is established and checked numerically, where 5 
governs the probability of return of a wave packet, pit) ~ , and D2 is the correlation 
dimension of the spectral measure at the critical point. 

Two applications which might be relevant for experiments are presented. First, we 
discuss the long time behavior of a wave packet of eigenstates with energies near E^. and 
show that in three dimension this behavior is drastically different from its two-dimensional 
counterpart Then, it is shown how rates for inelastic electron-electron and electron- 

phonon scattering are changed due to the critical dynamics as compared with the usual 
metallic case. We conclude with an overview of the results for various exponents in the 
two dimensional quantum Hall, the two dimensional symplectic, and the three dimensional 
orthogonal case. 



II. SCALING FORM OF D{Q,uj) 



The one-parameter scaling hypothesis ||12 



analysis of the metal-insulator transition |]13 



which underlies a great part of the theoretical 
allows to write the diffusion function D{q,uj) 



in d dimensions near the critical point in the form [14 



D{q,u;) 



2-d 



P 



-F{^/L^,^q) 



(6) 
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where ^ is the locahzation length, p = p{E) the density of states (DOS) near the critical 
energy E = Ec, -F is a scaling function, and 



:= iphoor'/'' = L{hujl^r^/\ (7) 

is a third length scale, besides and ^, relevant at the critical point (A = (L'^p)^^ is the 
mean level spacing). 

In addition to L^^, one can define a second frequency dependent length 



(8) 



the distance the particle diffuses in the time 1/uo. Here D{ijj) = limq^Q D{q,uj) is explicitly 
scale-dependent. Apparently, and scale with different powers of the frequency. How- 
ever, at the mobility edge this not the case. From the definition eq. (H), the dimensionless 
conductance for a hypercube of size L^^ is given by 

g* = G- = -a{u;)L^ . (9) 

Here, the frequency-dependent conductivity o-{uj) is related to the diffusion coefficient by 
the Einstein relation 

a{uj) = e^pD{uj). (10) 
Eliminating D{uj) from these equations leads to 

Thus, at the mobility edge both frequency-dependent length scales differ only by a factor of 

The general relation between the conductivity a and the conductance G for a sample 
in the form of a hypercube with sides of length L, G{L) = a{L)L'^~'^, and the vanishing of 
the /3-function dG/dL = at E^ fl^ requires the conductivity and therefore the diffusion 
function D ~ o" to scale like 

D ~ L^-'^ (12) 

at the critical point. 

A direct consequence of this scaling is a certain behavior of D{q, u) in different regions of 
the (g, ti;)-plane, as discussed by Chalker [|15|. In the limit of small q ^ 0, namely gL^^ <^ 1, 
sets the shortest length scale (region (A) in Fig. |l]), and 

D{q ^ 0, cu) ~ Ll-'' ~ cu^''-')/^ io-^O. (13) 

This behavior has been predicted by Wegner in 1976 [0. Recently, D(uj) ~ co^^^ has been 
confirmed numerically for d = 3 [|16i . 
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For larger q (region (B)), the length scale determining the scaling of -D(g, uo) is the inverse 
wave vector itself, if the localization length is not shorter than 1/q. This leads to 



D{q,uj) ^ q''-\ (14) 

It has turned out in the last years that there are corrections to the above scaling form 
due to strong fluctuations of the wave function amplitudes near the critical point. These 
fluctuations are observed numerically on small length scales corresponding to large g- vectors. 
In the (g, ti;)-plane, they give rise to another region (C) where D{q,uj), Eq. ([llD , acquires an 
additional factor g"'', the scale of which is given by either the localization length ^ or the 
length L^, namely 

D{q,u)r~.q''-'[q.mm{C,L^)]-\ (15) 
For our numerical analysis, we use the following scaling form of the diffusion function: 

gd-2 ^ 

D{q,u}) = -^f{x), X := [q-mm{^,L^)f , (16) 

which satisfies the scaling relation {C{E) > L^) 

D{q,uj) = b^-^D{bq,u}b'^) (17) 



which is required for a transformation of length scales r i— *■ r/b |jT4]. Note that in contrast 
to the conductance at Ec, the diffusion function itself is scale dependent. 

From these arguments it is not evident where the transition between regions (B) and 
(C) occurs. It is conceivable that instead of the scaling behavior in region (B) a smooth 
crossover from region (A) to (C) takes place. From the present numerical results we cannot 
address this question. 

The regions (A), (B), and (C) correspond to three different regimes of the function f{x): 
First, in order to have D{q —>■ 0,00) ~ L'^~'^ (regime A), the function f{x) must follow 
/(x) ~ jr{2-d)/d ^Y^e limit x 0. Second, for large x, one has / ~ ^"''/'^ in order 
to fulfill the behavior Eq. (|l^), regime C. The scaling form in the intermediate region (B), 
D[q,uj) ~ q'^~'^, requires that f{x) = const there. We thus can write the correlation function 
S{q, E, Lj) as 

S{q, E, uj) = x=[q- min(e, (18) 



77 



flu 1 + [xf{x)X 



where 



. ^{2-d)/d 

fix) = { 



X ^ (A) 

in between (B) (19) 
X — s> 00 (C) 



with constants Cq, cp and c^. In the brackets () we indicate the corresponding area in Fig. |I|. 
The boundaries between the three regions and the constants have to be extracted from 
numerical calculations. For the latter, we work close to E^. so that ^ > and the scaling 
variable x = [q- Ljf. The case C, < L^^ gives rise to the two remaining areas in Fig. |l] with 
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D ~ ^^"'^ (left to A) and D ~ (left to C), which are not relevant in our analysis. 

In the case of two dimensions, d = 2, regions A and B coincide and / is constant there. The 
numerical analysis (see below) shows that the large-x regime of f{x), governed by the power 
law involving the exponent rj, is indeed relevant over several decades of x. 

We notice that upon Fourier transforming S{q,uj) into real space, large values of x cor- 
respond to small r <^ min(^, L^^), where the strong fluctuations of the eigenstate amplitudes 
become important. In this region (corresponding to x = [mm{^, L^)/rY ^ 1), one finds 

S{r,u) ~ —^x-'+^Z'^ ~ ^[min(e,Lj/r]-'^+^ (20) 

which for cij ^ (,^ < L^) gives S{r, uj) ~ l/{uj^'^) ■ (r/^)^'', Eq. (|lD. For ^ > one recovers 
the form S{r^uj) ~ r~'^uj~^/'^ given by Wegner ||18|| . 



III. A WAVE PACKET AT THE CRITICAL POINT 

In order to illustrate the implications of critical dynamics, we use the relation between 
the diffusion function -D(g, oj) and the diffusion propagator P(g, uj), Eq. (^). The latter is the 
Fourier transform of the probability distribution P(r, t) describing the motion of an electron 
wave packet constructed from eigenstates with energies close to a fixed energy E. In the 
metallic region, P{r,t) is the solution of the diffusion equation in d dimensions 

^^M-W^P(r,t) = 0, t>0 (21) 
with the diffusion constant D, 

-ry{ADt) 

Due to causality, P{q, has to be analytical in the upper half plane as a function of complex 
uj: The integral 

1 f p-iuit 

P{q, t) = — duj . (23) 

has to be performed in the upper half plane for t < whence P(g, t < 0) = 0. On the other 
hand, for i > we can write 

P(q,t) = — I due-"^'\ ^- r^ + T r^l, t>Q (24) 

' 2tx J \-iu} + D{q,uj)q^ \-iuj + D{q,u)q^Y y ^ ' 

since for t > the integral is performed in the lower half plane where the second term gives 
zero contribution. For real tu, D{q,uj) is real and one can write 

1 r ,...~ic.t D{q,uj)q^ h I 



P{q,t> 0) = - / due-'-' '\ = - due-''-'S{q,E,u;). (25) 
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Recalling that S{q, E, u) is an even function in u and introducing the dimensionless variable 

X = {qL^Y = q'^ /ft'OJp, we obtain 

r°° r ^ / x^/ ^ N 2 f°° dx ( q'^t\ xfix) 
P(,,i>0)=2/^ ,.-cosM,E,.) = -l _ cos ^^1LL_. (26) 

Normalization of the wave packet requires P{q = 0,t > 0) = 1, i.e. 

2 f--dx xf{x) _^ ^27) 



vr Jo X 1 + 

which is an additional condition to be fulfilled by the function f{x) Eq. (|l^). Since the exact 
form of f{x) is not known, we cannot determine P{q,t) exactly. However, for large qH the 
probability P{q, t) in Eq. (|26D is determined by large x values due to the rapidly oscillating 
cosine-term. For large x, on the other hand, the form of f{x) is known, and we can write 



Here, we used the relation D2 = d — rj (see the following section). Furthermore, the lower 
(high frequency) cutoff x_ is determined by the inverse of a microscopic time scale r, below 
which the motion of the electron is ballistic and not described by the diffusion pole P{q,u). 

In the same line of argument, one obtains the probability of return of a wave packet to 
the origin in the long time limit. This quantity is defined by the r = value of P(r, t) and 
can therefore simply be recovered from P{q,t) by Fourier transformation. Again, the upper 
limit of the resulting g-integral is determined by a microscopic cutoff ~ 1// . One has 

Pit) = J^Jdq' nq. t) ~ dqq"-' {qH)-''-'' ~ t"^^/-^. (29) 

This dependence of p{t) is not valid for very short time scales corresponding to a (quasi)- 
ballistic motion of the wave packet, and times so long that the system boundaries become 
important. 

We can compare the critical behavior of p{t) in d = 2 and ci = 3 to the metallic case, 
where p{t) = {AnDt)-'^/'^: 



Pit) 




2d metallic 

2d critical , . 

3d metallic 
3d critical 



In particular, this demonstrates the different role of the exponent rj in two and three di- 
mensions: While in two dimensions, 77 essentially can be regarded as a more or less small 
correction to the ordinary metallic diffusion pole, the situation in three dimensions is drasti- 
cally different. There, due to the scaling form of D{q, u), the dynamical behavior of electrons 
near is qualitatively very different from the metallic case. 
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IV. RELATION BETWEEN THE CORRELATION DIMENSION D2 AND THE 

EXPONENT 77 



The relation Eq. (^) has been derived previously p]9| , |20| ; we shortly outline our proof 
since it comprises the definitions concerning the multifractal dimension D2 entering the 
analysis of the wave function at E^.. One first defines the probability 

P,(A) := / d''x\^E{x)\^ (31) 

to find the electron of fixed energy ii^ in a finite region (box i) of real space denoted by 
r2j(A) = as a subspace of the total volume VL = L'^, where / = XL, < A < 1. The scaling 
with A is used to define so-called fractal dimensions D^. One considers the averaged gth 
power of Pi(A) : 

Pi^, := ^ E [m)Y - A^"+('^-i)^^ (32) 

Here, A^(A) = A"'^ is the number of small cubes into which the original total cube of volume 
L'^ was split. Requiring that at Ec, the summation in Eq. (|32D over the different boxes 
labeled by i is equivalent to a disorder average in one (arbitrarily) fixed box (say i = 1), 



P(A,2) = (Pi(A)2), (33) 
one has (omitting the index i = 1 now) 

P{X,2)= [ d'^x I rfV (|^ij(x)H^ij(x')|'). (34) 

Because of the disorder average, the product of the squares should depend on the difference 
x — x' only. One can then introduce coordinates r = x — x' and R = {x + x')/2. The Jacobian 
of this transformation is 1, but because the integration region is finite, the integration limits 
are also changed which makes the evaluation tedious. However, for large enough Q{X) the 
integral over R simply gives the volume fl{X) itself, and we have 

P{X,2) = Q{X)- I d'r l\m Eim^^ E{r)?) . (35) 

Considering the definition of the two-particle correlation function Eq. (^ for u; — 0, one 
obtains 

P(X,2) = (Qp(E)?Q(X)- f d'^rS(r,E,uj ^0). (36) 

Jn(\) 

where Q is the total volume Q = Q{X = 1) of the system. Using Eq. (|l|), one has 

/ d^r-"^ ^ [^^ drr'^-^r-'^ (XL)'^-\ (37) 
Jnix) Jo 

Because Q{X) = {XlY, we find 
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P(A, 2) ~ n{X) ■ f d'^rr-" ~ A'^A^^''. (38) 

By comparison with Eq. ( p2D one reads off c? + (c? — 77) = Dq + D2- Noticing that Dq = d 
is the dimension of the total support of the wave function, we have the announced relation 
Eq. between the exponent 1] and the fractal dimension 

D2 = d-r]. (39) 

This fundamental relation relates properties of correlation functions to the spatial structure 
of the wave functions at the critical point. Furthermore, if rj is larger than at a critical 
point of a metal-insulator transition by Eq. (13) it follows D2 < d and the wave functions 



must be multifractal . 

We obtain another exponent relation by exploiting the long-time behavior p{t) of the re- 
turn probability of a wave packet, constructed from eigenstates near Ec. First, the exponent 
S introduced via p(t) ~ t^^ was generally proven by Ketzmerick et al. |^ to be equal to the 



generalized dimension D2 of the spectral measure. Just as D2 describes spatial correlations, 
D2 describes the correlations of the local density of states of the system as a function of 
energy. In |^ it was shown that 5 = D2- On the other hand, we have seen that 5 = D2/d, 
Eq. (^91), and therefore 

S = £l=D2. (40) 



V. NUMERICAL ANALYSIS 

A. Model for the numerical investigation 

The dynamics of non interacting electrons in the presence of disorder is studied within 
the framework of the Anderson model described by the Hamiltonian 

n = Y.e,\v){v\ + K-,r'|r)(r'|, (41) 

r (r,r') 

where the diagonal disorder potentials er are independent random numbers with a constant 
probability distribution in the range —W/2 < < W/2 and the non diagonal transfer 
matrix elements between nearest neighbors, Vry, are taken to be the unit of energy. The 
vectors r denote the sites of a simple cubic lattice with lattice constant a and periodic 
boundary conditions are applied in all directions. 

Eigenvalues and eigenvectors have been obtained for systems of size up to {L/a = 40)^ 
sites by direct diagonalization using a Lanczos algorithm. The critical regime in the middle 
of the disorder broadened tight binding band, where both the localization length and the 
correlation length diverge, is known to correspond to a critical disorder of Wc/V ~ 16.4 
which separates localized {W > Wc) from metallic behavior {W < Wc). 
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B. The function Z{uj) and the exponent fi = rj/d 



For the numerical analysis, we defined the function 

Z{E,E'):= J d'^x\^E{x)\'^\'^E'{x)\'^ (42) 

for eigenstates with energy E and E', where d denotes the spatial dimension. In the following, 
we will be interested in the case where the localization length is the largest length scale in 
the system so that C,{E) > L^. We show that one can directly extract the quantity r] from 
Z{E,E') which is easier to obtain numerically than, e.g., a direct determination from the 
diffusion function D{q,uj) itself. Indeed, the two-particle correlation function in real space 
for r = 0, 

S{r = 0,E,uj) := (j26{E- - E^)5{E+ - E^)^s-(0)^^,-(0)^s+(0)^^+(0)\ (43) 

\a,/3 / 



can be related to Z[E, E') assuming that the disorder average in Eq. ( P3| ) is equivalent to a 
spatial average for one fixed impurity configuration. In this case. 

Sir = 0,E,u;) = {np{E-)) (^^p(^+)) J d'^x \^ eAx)\'\'^ e+{x)\\ (44) 

with huj = E^ — E~ and E = [E^ + E^)/2. The equivalence of these averages is underlying 
the subsequent analysis of the function Z{E,E'). The relation 

obtained by Fourier transformation, is used in the following. Keeping E fixed at the critical 
energy, the quantity Z depends on u = {E^ — E~)/h only. Since the numerical calculation 
is performed on a finite lattice, the g-integrals have to be cut off at g = 27r/a, where a is the 
lattice constant. A lower cutoff is given by the system size L itself. In terms of the scaling 
variable x, the lower cutoff is x^{u) := {{L/2nyphu)~^ = {2'7iYA/huj. The upper cutoff is 
x+{u!) := {{a / 2t[Y phuj)~^ = X-{uj){L/aY. In a finite system, energy transfers huj from one 
state with energy E to another with energy E + huj satisfy hu ^ A. In terms of the variable 
X, this means x^{uj) ^ (Svr)*^. For the function Z{uj), Uuj = E^ — E~ in Eq. (^), and we 
obtain 

Z(u}) = , f , r''^ dq q'^Siq, E, u) = , H^^^ dx ^/^f,^ (46) 

Here, 5^=2 = 27r and 5^=3 = 47r. For small uj and large L the integral is dominated by the 
large x behavior of f{x) (region (C)). 



For small cj, there is a range where Z{uj) ~ uo^'^ with fi = T]/d. A logarithmic plot of 
Z{uj) thus yields the exponent /i = T]/d, furthermore from Eq. (^) one can determine the 



coefficient c^. 
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1. Numerical data 

We used the numerical data to obtain the exponent fi = rj/d from the function Z[uj). 
The result for a system of size Vl = (40 a)^ is shown in Fig. ^. One clearly observes a power 
law with an exponent fi = 0.5 ± 0.1. This yields rj = fi ■ d ^ 1.5. In the energy range 

^ 300 A > |E - E'l > A, we can fit the function Z hy Z{uj) = Zo\E - E'\-°-^ with 
Zq = 2.5 ■ 10~^. A/V is about 2.7 ■ 10"'^ in the present system. ' For even larger values of 
uj or correspondingly smaller x, one should enter regions (A) and (B) which are beyond the 
applicability of eq. (|47|). 

On the other hand, for small u ^ A, Ziuj) saturates at the value of the inverse partici- 
pation ratio P^'^\E) = J d'^x\^E{x)\^ oc (L/a)^^^. A similar limiting behavior has also been 
observed in the QHE-case (see below) and in the 2d-symplectic situation 

2. Comparison to the quantum Hall case (d = 2) 

In the quantum Hall case there exists no complete Anderson transition, but for finite 
systems a critical behavior and multifractal eigenstates [0-0 can be observed in an energy 
range about the center of the disorder broadened Landau band, Eq, where the localization 
length exceeds the system size, ^ = ^q\E — Eq\~^ » L. We used our previously obtained 
eigenvalues and eigenfunctions of a system of size 125 x 125 and a magnetic field which 
corresponds to 1/5 flux quanta per plaquette (see for details) to calculate the function 
Z{uj) shown in Fig. ^. Again, a power law relation can be observed with a saturation at 
UJ ^ A = 3 ■ 10^^ V. The exponent fi = r]/2 = 0.26 ± 0.05 is somewhat larger than the value 
obtained from the exponent 6 of the temporal decay of a wavepacket {ri/2 = 1 — 6 = 0.19) 
built from the same critical eigenvectors [Q. 

C. Exponents D2, D2, 6. 

1. The Wave packet 

The correlation function C(t) related to the probability of return p{t) |2^J§] of a wave 
packet, constructed from eigenstates with energies E ^ Ec of a 3d Anderson model, is shown 
in Fig. ^ where C(t) is plotted as function of tA/h. A power law p(t) ~ can be clearly 
identified, we obtain S = 0.6 ± 0.05 in the region where t > 2 ■ lO^^h/A. The condition 
qH/fip in Eq. ( p8l) can be checked if we use q = 2tt/L as the smallest q value in the system 
{L/a = 35). Together with the density of states, which is p = 0.058/(V^a'^), we obtain from 
Eq. (|28|) the condition t 3> 4 ■ 10~^h/A which is in accordance with the observed behavior 
in Fig. |. 

2. Multifractality exponents 

We calculated the correlation dimension of the multifractal wavefunction D2, Eq (P^, 
and the correlation dimension of the spectral measure, D2, defined by 
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7(g,e) = liml/L2EE( |a^p)''~.(-l)^^ (48) 

using a box-counting method. The result of the scahng behavior of the latter is shown in 
Fig. 1^. The exponent 6 from the probability of return of a wave packet (Fig. ^) and the most 
probable scaling exponent, ao, determining the characteristic /(a)-distribution (see e.g. ||2^ 



and references therein) related to the generalized fractal dimensions Dg which completely 
describe all the moments of the spatial amplitude fluctuations of the critical wave functions 
have also been calculated. Our results are compiled in Table 1. The relations concerning 
the exponent rj which have to be fulfilled according to Eq. (^), Eq. (^), and the definition 
/i = r]/d, are r] = d{l - 6) = 1.2 ± 0.15, r] = d - D2 = 1.3 ± 0.2, and r/ = 1.5 ± 0.3, 
respectively. Thus, one can say that within the numerical uncertainty the different methods 
of determining rj yield the same value 1] ^ 1.3. We note that our value for D2 is very close 
to a result of previous numerical work by Soukoulis and Economou using a different method 
|1|] who obtained D2 = 1.7 ± 0.3. 

VI. INELASTIC SCATTERING AT THE ANDERSON-TRANSITION 

In this section we investigate the implications of critical dynamics at the Anderson tran- 
sition on inelastic scattering rates. Inelastic scattering rates were calculated previously ||24 



for a quantum Hall system where the critical energy coincides with the Landau band center. 
It could be shown that the exponent r] describing the eigenfunction correlation showed up in 
the temperature dependence of the scattering rate and the energy loss rate of electrons with 
acoustical phonons. On the other hand, the temperature dependence of the electron-electron 
scattering rate was similar to that in a two-dimensional disordered metal without a mag- 
netic field. It is therefore of some interest to study these quantities for the three-dimensional 
system, too. 

A. Electron-phonon scattering rate 



The electron-phonon (e-p) scattering rate is calculated in the standard way [31| to second 
order in Vg, the e-p coupling matrix element. It is assumed that the impurity potential is 
not changed by the lattice motion. In a metal this assumption is not fulfilled a priori: The 
impurities are embedded in the lattice and therefore move in phase with the other lattice 
atoms |^2[. Here, we address to a situation where the random potential is generated by 
external sources like donor atoms far away from the electron gas. 

The imaginary part of the self energy gives the inelastic lifetime of an electron of energy 

E 

rr,\E) = ^Y.\Vf\MEM^)\'^ (49) 

I^E' 

X {[n{ujg) + f{E')]6{E -E' + hUg) - {Ug ^ -Ug)} 

Here, n{ujg) denotes the Bose function for phonon frequency Ug and / the Fermi function. 
The terms in the curly brackets correspond to phonon emission {+hujg) and -absorption 
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{—hujq), respectively. The £"-sum runs over all intermediate eigenstates of energy E' which 
appear in the momentum matrix element Me^e'{^) = {E\ exp(iqx)|£"). The central quantity 
containing the information about the unperturbed system in ( ^Q]) is 



\MEM^)\^SiE -E' + hujg) = Qp{E + hujg)\ME,E+n.M)\^, (50) 

E' 

where p is the three-dimensional density of states and Q is the total volume of the system. 
We express the matrix elements by the two-particle correlation function 

np{E + nuj)\ME,E+U^)\' = -^fi^-^ E + ^,u) (51) 

We restrict the discussion to energies E = Ee, where Ee is the Fermi energy. Using 
n{ujq) + f{E + huiq) = 1/ sinh(/3^a;q), we get 



^'P~27chJo "^"^ sinh Phco "''p(i?)" ' ^^^^ 

where we defined the function a'^F{u!) = o?iJ^ for acoustical phonons with coupling constant 
o?. In the case of non-piezoelectric materials, only the deformation potential coupling is 
relevant. Then, n = 3 and = (E'^h)/ (2pMC^s) with speed of sound cs, deformation 
potential S, and mass density pm- The temperature dependence at the critical energy 
E = Ec is determined by the form of the correlation function, defined via Eq. (p!9| ) with 
d = 3. The phonon dispersion uj = csQ leads to a; = u"^ /{c^gph). The dimensionless variable 
X ((i = 3) can be written as x = {uj/ujcr)'^ by introducing the crossover frequency 



{hclpf\ (53) 



The exponent for the T-dependence of r^^ can be extracted in two limits: For frequencies 



UJ ^ ojcr one has x ^ 1 and the correlation function S ~ {xf{x))^^/uj ~ c<j^^+^^/^. In 
the opposite case uj -C tUcr and x <^ 1, S* ~ (x/(x))/co' ~ uj^^^~^ . This leads to different 
behaviors of the scattering rate T~p in the limit fc^T <C hujcr and ksT ^ Tiujcr-, respectively. 
We can evaluate Eq. (|52D, introducing (3?iuj as new variable. First, for temperatures fc^T ^ 
huE) ^ hucr larger than the Debye energy, the dependence is linear in T. This is (as in an 
ordinary metal) the trivial high-temperature case. The more interesting limit ksT <^ Huje 
gives, together with a'^F{uj) ~ cu" and the knowledge of the correlation function S, the 
limiting forms 

^ev ~ I 3^n-2+2^/3 . ^^J. ^ ^^^^^ 1^34 j 

These two regimes correspond to the case (A) and (C) in Eq. (|T9|). The full temperature 
dependence of can be obtained from a numerical evaluation of Eq. (^). 

The result r^^ ~ 2^3+4/3 f^j, deformation potential scattering [n = 3) in the low tem- 
perature limit is a direct consequence of the 'Wegner-scaling' of the diffusion function 
D(q,iu) ~ uj^^^ in regime (A). It should be compared to the dependence r'^^ ~ for 
the corresponding rate in the disordered metallic case We also note that as in the 



quantum Hall case [Q, the exponent r] appears only in an intermediate and not in the 



low-temperature regime. 
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B. Electron-electron scattering 



Electron-electron (e-e) scattering rates at the Anderson transition have been calculated 
first by Belitz and Wysokinski ||3^ , who used the so-called 'exact eigenstate formalism' and 



an RPA-like approximation for the Coulomb interaction. The advantage of this method is 
that it is non-perturbative in the disorder and, at least in principle, applicable to a large 
number of systems. By this one can generalize results for e-e scattering rates which are 
already known in the ballistic or the diffusive, weakly localized case. However, one has to 
keep in mind that this approach is perturbative in the interaction between the electrons. 
The general expression for the electron scattering rate at finite temperature ksT = 



1/(5 > 0, given by the on-shell electronic self-energy at Fermi energy E = Ep ||3^j33 

0(^) 



r7^ = 2 du 



e 



oo 



sinh jShuj 



'^(^) = ir?T Im V(g, ^)Im <l>(g, uj) (55) 

nnp\l q 



can be derived by starting from the exact eigenstates of the unperturbed system Here, 
uj) is the density relaxation function of which the imaginary part is related to that of the 
density response function u;) by Im$(g, cj) = (l/a;)Imx(g, cj) = nS{q] E,u) where the 
energy E appears in the density of states p = p{E). Notice that Im$(g, tu) = Im$(g, —a;) 
and lmV{q,uj) = —lmV{q, —u) because Rex{q,^^) is an even and Imx(g, tu) is an odd func- 
tion of u!. One therefore can restrict the integration to positive values of uj. The above 
formula can be viewed as the rate for the scattering of one electron at the fluctuations of 
the electromagnetic field caused by the motion of all the other electrons the function 
being an effective density of states of these fluctuations. It corresponds in the case of 
electron-phonon scattering to the Eliashberg-function a'^F which there essentially gives the 
phonon density of states. In the Coulomb potential, screening is included via the RPA 
approximation 

J'^"] (56) 

with the bare Coulomb potential VQ{q) = 47re^/g^. Notice that only the inclusion of dynam- 
ical screening leads to an non vanishing imaginary part of V{q, u) and thus to a nonzero e-e 
scattering rate. The bare Coulomb potential would render the electron self-energy merely 
real and t~^^ would be zero. The influence of arbitrary disorder on the screening here is 
incorporated in the correlator x which contains all the information about the unperturbed 
system. In the g = limit of the diffusion function D{q = 0,u)) was used for the 
evaluation of Eq. (|55| ) at the critical point. Although the integration requires both the q- 
and the w-dependence of the diffusion function, it turned out that this approximation affects 
only the prefactor of the scattering rate and not its temperature dependence. 

In the following, we will use the full form of the diffusion function as it enters into the 
density-density correlation function, which at the critical point can be written as 

X f (x^ q'^ 
X{q,u;) = ip-———, x = - — , (57) 
l + zxj{x) hup 
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The function ^{u) can than be evaluated; its zero- frequency hmit is given by 



$(cj = 0) = $0 



dx 



6tt^ Jo 1 + {xf{x)Y 



(58) 



As in the two-dimensional diffusive and Quantum Hall case |36|,^ , the dimensionless quan- 
tity $0 can be expressed by an integral over a function of the variable x = j {hup) where 
d is the dimension. As in the above cases, the scattering rate according to Eq. (|55|) formally 
diverges because the function $(a;) describing the DOS of the electromagnetic fluctuations 
from the electron 'bath' does not vanish at zero frequency. However, for finite temperatures 
such that the scattering rate is larger than the frequency uj, the description of the critical 
dynamics in terms of xi'f 1^) can no longer be valid and the integral should be cut off at 



31 . This leads to a self-consistent equation (for low T, $(u;) can be replaced by $o 



r. 



-1 



2$ 



oo ]^ 

duj . , ^. , 

-1 smh pnuj 



of which the solution can be written as 



h 



-7, 



7 



-2<l>o lntanh(7/2). 



(59) 



(60) 



This is in fact the result of Belitz and Wysokinski, differing only in the pref actor given 
by 



V3 

7 = lntanh(7/2) 

4c 



Ref. [33], 



(61) 



where c = 13.4 or 1 [^], respectively. In view of our analysis, Eq. (|T9|) and Eq. (^8l), we 
find that the prefactor depends on the microscopic details of the system, namely the con- 
stants Ca, c^, and the dimensionless conductance at the critical point g*. In two-dimensional 



quantum Hall systems the critical conductance was found to be universal |T^. However, we 
know of no argument that shows the critical conductance to be universal in three dimen- 
sions. Thus we find that the universality of the rate claimed in is related to the 
universality or lack thereof of the critical conductance g* . 



VII. CONCLUSION 

We have discussed some aspects of multifractal fluctuations near the Anderson transi- 
tion in three-dimensional electron systems. The behavior of the dynamical diffusion function 
D{q,u!) was reviewed. For large wavevectors q and small frequencies u multifractal fluctua- 
tions lead to the occurrence of an anomalous diffusion exponent t] in the diffusion function. 
We have found ^ 1.3 in a numerical analysis of eigenfunctions obtained by numerical diag- 
onalization, and verified that t] is related to the correlation dimension D2 of the multifractal 
fluctuations hj i] = d — D2. 

The multifractal eigenfunction fluctuations influence the temperature dependence of the 
electron-phonon and electron-electron scattering rates. We found that the low-temperature 
exponent of the former is modified by the finite value of rj, while the latter is linear in T 
and depends on the dimensionless conductance. 
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TABLES 



System 


V 


b 


"0 




^2 


11 = rj/d 


2d QHE 


2.35 


0.81 ± .02 


2.29 ± .02 


1.62 ± .02 


0.80 ± .05 


0.26 ± .05 


2d Sympl. 


2.75 


0.83 ± .03 


2.19 ± .03 


1.66 ± .05 


0.83 ± .03 


0.175 ± .03 


3d Ortho. 


1.45 


0.6 ± .05 


~ 4.0 


1.7 ± .2 


0.55 ± .05 


0.5 ± .1 



TABLE L The calculated exponents related to the critical states at the metal-insulator transi- 
tions in two and three dimensional systems. The results for the critical exponent of the localisation 



length v are taken from |Q (QHE), (2d symplectic), and [26,2^ (3d orthogonal). The re- 
maining values for 2d systems not calculated in the present work were taken from (QHE) and 
(2d symplectic). 
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FIGURES 




FIG. 1. Five different regimes (separated by solid lines) for the diffusion function D(q,uj) 
after Chalker [|l5|. Relevant for our numerical analysis are the three regimes A,B,C — > oo) 
with the behavior of D{q,uj) and in particular q'^ "^{qL^) ^ with the exponent 

r/ = D2 — d related to the multifractality of the eigenstates. According to the value of the scaling 
variable x = {qLi^Y, the function /(x), Eq. (|l9|), follows a different power-law in x. The dotted 
lines indicate this crossover between A, B and C. Since q^ = (a) corresponds to small, 

(b) to intermediate and (c) to large x. The latter case is relevant for the numerical determination 

of T]. 
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FIG. 2. Energy correlation function Z{io) of eigenstates taken from tlie critical regime 
{Wc/V = 16.3) of a 3d Anderson model of system size {L/a = 40)^. We used 371 eigenstates 
with energies E from the interval —0.05 < E/V < 0.05 to extract the exponent // = rj/S = 0.5 
from the power law decay. 
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FIG. 3. Energy correlation function Z{uj) for a QHE-system of size {L/a = 125)2. 330 eigen- 
states taken from an energy interval /S.E = Q.IV around the center of the lowest Landau band 
were used to determine the exponent jj, = rj/2 = 0.26. 
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FIG. 4. Temporal autocorrelation function C{t) for a 3d critical Anderson model 
iWc/V = 16.5) of size {L/a = 35)^ versus time t in units of the Heisenberg time tn = h/A. 
To construct the wave packet at the initial site a total of 1194 eigenstates from the energy interval 
[—0.5,0.0] have been used. 
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FIG. 5. The correlation dimension D2 = 0.55 it 0.5 of the spectral measure as obtained from 



the scaling relation of 72 (see Eq. 48) versus the length of the energy interval e for a critical system 
of size {L/a = 40f. 
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